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I. Introduction ; 


Several Investigations wore carried out. In the following I 
shall give: a list of publications, highlights of the findings, 

and research in progress not yet written up for publication. 

2. Pub I i cat ions : 

The rese< rch resulted in the following papers which have been 
submitted for publication: 

a: Robinson, W.O., and R.C. Srivastava, 1981: Calculations 

of Hailstone Growth in a Sloping Steady Updraft. Submitted 
to Atmosphere - Ocean. Also LAP TR . 

b: Srivastava, R.C., 1981: A Simple Model of Particle Coalescence 

and Breakup. Submitted to the Journal of Atmospheric Sciences. 
Also LAP TR . 

c: Matejka, T.J. and R.C. Srivastava, 1981: Doppler Radar Study 

of a Region of Widespread Precipitation Trailing a Mid-Latitude 
Squall Line. To appear in the Preprints of the 20th Conference 
on Radar Meteorology. (A more detailed paper will also be 
submitted for Journal publication.) Also LAP TN #10* 

It may be mentioned t^at most of the work for (a) was done under an 
NSF grant. The NASA grant provided the opportunity and support for organizing 
the results of the research and writing it up for journal publication. The 
work for (b) was supported by NASA and NSF. The data collection on which (c) 
is based was supported by NSF; NASA and NSF provided support for the data 
analysis and research. 



3- Resul ts 


Full details of the results of the investigations may be found in 
the above papers (copies attached). In the following I shall sunvnarize 
the main results and their signi f icince. 

• Hall Growth ; 

Sustained falls of large hail are usually associated with clouds 
having a sloping updraft which can be regarded as quasi^steady . A 
number of studies have demonstrated that the trajectory of the growing 
hall is important in considerations of hall growth and artificial modi- 
fication of hall processes. 

A simple model of a sloping updraft and hail growth was set up for 
which analytical solutions could be obtained. It was found that the 
.starting position of the hall embryo has an important effect on its 
subsequent growth. Three regions of embryo starting positions were 
distinguished,. Embryos starting in the central portions of the updraft 
are carried up rapidly and displaced far from the updraft into the 
anvil. Most of these particles cannot survive during fall to the ground. 
Embryos starting near the upshear edge of the updraft have a simple 
arcing up and down trajectory and produce hailstones most of which can 
be expected to melt before reaching the ground. Embryos starting in 
the third region, namely, the region near the downshear edge of the 
updraft give rise to the largest hailstones and have predominantly 
re-ci rculating trajectories. 

It is believed that the subject deserves further study; an investi- 
gation of the general properties of particle trajectories in model flow 
fields can give insights into hail growth processes which will be diffi- 
cult, if not impossible, to obtain from numerical models. 


3.2 Particle Coalescence and Breakup 

The processes of particle coalescence and breakup^ spontaneous and 
collisional, are important in determining the evolution of particle 
size spectra. A firm understanding of the evolution of particle sizes 
is important for a number of reasons: The particle size distribution 

is involved in various nwthods for measuring precipitation rates and 
in numerical studies of cloud dynamics. 

Previous studies have identified, to some extent, the roles of 
the various processes in shaping particle size spectra. Coalescence 
tends to produce particles of larger sizes and exponential spectra 
of the kind observed In nature. Coalescence and breakup tend to 
produce stationary exponential distributions. The effects of 
collisional disintegration dominate over those of spontaneous disin- 
tegration as the precipitation content increases, and the equilibrium 
spectra tend to be parallel to each other. 

The above findings from numerical models are verified by analytical 
solutions of the kinetic equation for the particle size distribution 
for simplified forms of the collection and breakup functions in paper (b). 
The exact analytical solutions provide insights into the processes 
governing the forms of particle sire spectra not available from numeri- 
cal calculations. 

3.3 Precipitation Associated with a Squall Line 

During PROJECT NIMROD, extensive Doppler radar observations of squall 
lines were obtained in Illinois. Paper (c) discusses observations of a 
region of widespread precipitation trailing a squall line. Such regions 
of widespread precipitation are often observed in association with mid- 
latitude squall lines and more generally in the case of tropical squall 
lines. Physical considerations suggest that the region of widespread 
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precipitation may have meso-scale vertical motions In It which could be 
important dynamically as well as for budgets of various quantities. An 
important budget quantity is the rainfall. It is often observed that 
while the intense rain associated with the convective cores of a squall 
line is relatively spotty and short-lived, the light rain associated 
with the region of wide-spread cloudiness last for a long time and extends 
over a much larger area. It has been estimated that the precipitation 
associated with the region of widespread cloudiness often contributes a 
significant fraction (=0.5) of Che total precipitation from the squall 
line system. This also supports the idea that vertical motions exist 
In the region of widespread precipitation. 

In paper (c) , we have extended the VAO method to obtain accurate 
estimates of small vertical motions (w = 10 cm s"^) and applied it to 
observations of widespread precipitation associated with a squall line 
system. The results show generally ascending motions in the upper parts 
of the cloud and descending motions In the lower parts of the cloud. 

It is suggested that precipitation evaporation associated with the 
descent, and precipitation production associated with the ascent help 
maintain the observed verticas air motions. 

4. Research in Progress 

Work discussed in 3.3 is continuing. We anticipate obtaini.ig 
quantitative estimates of the rain output from the squall line system 
using data from a dense network of raingauges (average spacing = 3 km). 

5. Copies of the papers mentioned in Section 2 are attached. 
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Abstract 


Analytical solutions to the equations of growth and motion of hailstones 
in updrafts and cloud water contents which vary linearly with height are 
used to investigate hail growth in a model cloud. The model storm is 
steady and two'dimens ional and is constructed by compositing sections with 
linear variations of vertical air velocity and cloud water content so as 
to approximate the Fleming storm. Hall embryos are introduced at various 
positions in the updraft and their subsequent history calculated. A 
strong correlation is found between the embryo starting position and its 
trajectory and final size. Embryos starting in the central portions of 
the updraft are carried up rapidly and displaced far from the updraft. 

Embryos starting near the upshear edge of the updraft have a simple arcing 
up and down trajectory and produce small hailstones which can be expected 
to melt before reaching the ground. Embryos starting near the downshear 
edge of the' updraft give rise to the largest hailstones and have predominant! 
recirculating trajectories. Effects of changing some of the model parameters 
are investigated. 


INTRODUCTION 


Th« foundations of the theory of ha 11 growth were laid by Schumann 
(1938) and Ludlam (1958). Ludlam showed that the requirement of heat- 
balance does not always allow all the water accreted by a hailstone to 
freeze. In Ludlam's model, the excess (unfrozen) water was shed; this 
limited the maximum size to whicii a hailstone could grow to a diameter 
of about 4 cm. However, experiments by List (1959) showed that the unfrozen 
water need not be shed from the growing hailstone but rather that the hail- 
stone can grow in a "spongy" mode. By assuming that no water is shed from 
a growing hailstone, Hitschfeld and Douglas (1963) showed that large hail- 
stones can grow in a one-dimensional (vertical) cloud, 'n times consistent 
with observations, provided high concentrations of water are present as 
postulated by Marshall (I96I). In their work on hail growth, Sulakvelidze 
et a 1 . (1967) also considered regions of high water content or accumulation 
zones , 

Building on the above pioneering work, and spurred by possible appli- 
cation to artificial modification of hail, a considerable volume of research 
on hail growth has been reported in the recent literature. The simplest 
models of hail growth are one- di mens ional models in which the vertical air 
velocity and cloud microphysical properties are prescribed functions of 
height or height and time. (e.g. List et a 1 . , I968; Mus i 1 , 1970; and Dennis 
and Musil, 1973). A somewhat more complex but still a one-dimensional model 
of hail growth considers interacting dynamics and microphysics (e.g., 


Wisner et a 1 . , 1972 and Danielsen et al . , 1972). 



A criticism of one-dimensional models is that the storm dynamics are 
not adequately modeled or simulated in one dimension. Further, observations 
show that, at least in middle latitudes, storms producing sustained large 
hail occur in an environment having a vertical shear of the horizontal wind, 
suggesting that the horizontal dimension plays an important role in the 
hail process. This has been stressed by Ludlam (19^3) and others, 

Calculations of hail growth in two-dimensional nwdels were reported 
by English (1973), Orville and Kopp (1977) and others. Guided by certain 
observed characteristics of Alberta hailstorms, English prescribed the 
vertical air velocity and water content as functions of vertical and hori- 
zontal distance, and showed that hail of observed sizes could be grown in 
the clouds in reasonable times. The liquid fraction of the hailstones 
grown in the nradel were reasonable and not excessive as in the one-dimensional 
model of Hitschfeld and Douglas. Orville and Kopp considered a two- 
dimensional mode! with highly parameterized microphysics. It is doubtful 
that the dynamics are properly simulated in two dimensions and the calcu- 
lation of hail growth suffers because of the extensive parameterization of 
microphysic.". However, the simplifications of two-dimensionality and the 
microphysical parameter i zat ions were dictated by computer limitations. 

It appears that the dynamics, microphysics and the horizontal and 
vertical dimensions all are important in the problem of sustained falls of 
large hail. A detailed numerical calculation considering all the parameters 
cannot be performed on present-day computers. Also it is doubtful that such 
a calculation would be very illuminating. Therefore, qualitative models 
based on physical considerations ar.d observations become important. Among 
such models, we mention those of Browning (1963, I96A, 1977) and Browning 




and Foote (1976) which clearly demonstrate that both the vortical and the 
horizontal dimensions are important in tha process of hall growth. 

In this paper, we present a two-dimensional, steady-state model of 
hail growth In a sloping updraft which Is simplified to the point that 
analytical solutions for hailstone growth and traJec*’ory become possible, 

In spite of its simplicity, the model 1 1 Iul ♦■rates the effects of particle 
trajectory on its growth and reproduces certain aspects of the qualitative 
models mentioned In the previous paragraph. The model used here is an 
extension of the two-dimensional model of precipitation growth in an 
upright axial ly-symmetric steady-state updraft presented by Srjvastava 
and Atlas (1969). Other two- and three-dimensional models which compute hail or 
precipitation particU^ trajectory and growth are those of English (1973), 

Sartor and Cannon (1976) and Paluch (1978), 

2. THE MODEL 

We consider a steady-state two-dimensional model of a storm. The 
model storm consists of a number of "sections" in the horizontal and the 
vertical in each of which the vertical air velocity and the "effective" 
cloud water content are linear or piece-wise linear functions of the height 
and constant in the horizontal direction. Discontinuous changes in certain 
parameters are permitted between horizontal sections. These assumptions 
allow analytical solutions to be obtained for the particle size and position. 

We shall now present the equations for the model. 

2. I Equations 

The following equations apply to a portion of the model cloud in which 


the cioud properties vary linearly with height and are constant in the 




horizontal. Details of some of the derivations can be found In Srivastava 
and Atlas ( 1^0:^) . 

The vertical air velocity Is given by 

w » Wq + az ( 1 ) 

where and a are constants and z Is the vertical coordinate. The "effective" 
liquid water content for hail growth is considered tc be the product of 
M, the cloud water content, and E, the collection efficiency of hail; 

Wf? assume 


ME - Yq + Y jZ (2) 

where Yq and Yj are constants. The terminal velocity V of hail of diameter 
0 is taken as 


= KD (3) 

where K (= 2 x 10 cm s is a constant. This should be acceptable for 
i mm. The continuous collection equation for hail growth is 


^ _ ME 
dt " 2p 


ik) 


where p (taken f.s 0.9 9 cm ^) , is the density of the hailstone. It is 
assumed the t no water is shed from the hailstone, A weakness of this 
(nodel is that the liquid fraction of the hailstone is not calculated. 
The vertical position of the particle is given by: 


dz 

dt 


w V 


(5) 


The above equations can be combined to yield the following equation for 
the z and V of a growing particle: 





z = -p 


o 


d^z 

dt^ 


( 6 ) 



where 


are constancs. Analytical solutions for the above second order linear 
differential equations with constant coefficients can be easily obtained. 

We now consider the horizontal Position x of the particle. It is 
given by 


where u is the horizontal velocity, 
. 3 It 


3 (p w) 

a 


'a 3z 


3 z 


For slab symmetry u is given by 

( 10 ) 


The air density p is assumed to be 
3 

p » p exp(-Az) (1 1 ) 

a O 

where A is a constant, we take 

A - ( ^ - f)/T (12) 

where g is the acceleration due to gravity, R is the gas constant For air 
and Y the mean lapse rate of temperature and T the mean temperature. 

From the above equations, we have 


dx 

dt 




X- [Xw 


a] 


(13) 


The above equation can be Integrated to give the trajectory of an air 
particle starting at (x^, z^) 


w +az, 

= X • — 2, .exp[X(z-z )] 

o w +az ^ o 


X 


(14) 



and the position of a growing hail particle 

X 3 a , 

X * X exp[(?tw a - a )t + I'V-V )] , 3 , 0 

o o P j P J O I 

A3 

^o “ ~ ° (’5) 

where V is the initial fall speed of the particle at (x , 2 ) . 

The solution of ( 6 ), (7) and (l*>) gives the particle fall speed and 
vertical and horizontal positions in terms of the initial conditions 
(V^, z^) and the nxjdel parameters (w^, a, 3 ^, 3 p X) . 

2.2 Model Composition 

The cloud described in the equations of Section 2.1 does not have 
enough flexibility to n»del realistic sloping updrafts. The airflow to 
the right of the plane of symmetry of such a model is shown schematically 
in Fig. 1 . At the lower levels, air converges inward and at the higher 
levels, the air diverges outward into the anvil. By appropriately 
piecing together portions of the model of fig. 1 , with other portions of 
this model, or portions of models with different parameters, it is possible 
to approximate fairly general sloping updrafts. Thus we can take the portion 
ABCD of fig. I and join It with a portion from a different model at the 
top or the sides, and continue the process to get a composite model. Cer- 
tain constraints must be observed in this compositing process. Continuity 
of mass must be considered at the interfaces. The streamlines (Eq. lA) are 
such that no air crosses them. Thus the choice of left and right edges of 
a section is the streamline at the edge. It can be seen from Eqn, (lA) 

that the streamlines do not change i f w /cx remains constant. Thus sections 

o 
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Figure 1. 

Sohematio air 'particle trajeatoTriea 
in a slab ayrmetria aloud. 



Figure 2. (a) The. updraft shape, and the distributions of (b) vertical 

air velocity and (c) effective alo^ld water content in the 
model Fleming storm. 




can be placed side by side if they haiVe the same value of w^/a. Still air 
sections have no streamline constraints. !n piecing sections vertically, 
we require continuity of mass flux. The vertical mass flux through a 
horizontal boundary is wL, where L is the length of the horizontal boundary. 
In this model, vertically adjacent sections have the same L; thu*: the verti- 
cal velocity is required to be continuous at the horizontal interfaces 
between sections. 

2. 3 Mode I pa; ameters 

The parameters of the model storm will be selected to approximate 
the quas i -steady-state phase of the Flemin9 Storm (Browning and Foote, 

1976). The model storm consists of three adjacent sections (Fig. 2), 

The updraft velocity in the three sections is shown by the correspondingly 
numbered curves in Fig. 2b. Section 3 is composed of still air. Section 
2 is a transition region between the core of the updraft^ Section I, and 
the still air section. The vertical velocity in Section 2 is 1/3 of that 
in Section 1, In each of the sections, the updraft and effective liqu'd 
water content are piecewise linear functions of the vertical coordinate and 
constant in the horizontal. The vertical air velocity in the core of the 
updraft near the cloud base is 11.2 m s ^ and is coinparable to the vertical 
air velocity measured by instrumented aircraft near the base of the Fleming 
cloud. The vertical velocity in the model is zero near the observed top of 
the Fleming cloud. No measurements of the vertical air velocity were 
available in the body of the Fleming cloud, Using parcel theory (with no 
water loading) the temperature sounding data yielded the updraft shown by 
the continuous curve. The model updraft assumes a lower maximum vertical 
air velocity located near the middle of the cloud rather than near the 



cloud top. According to Sulakvelide et al . (196?) and Marwitz (1972), 
the maximum vertical air velocity occurs in the middle levels of cloud. 

The assumed effective liquid water content in the model cloud 
(fig. 2c) is an approximation to the computed adiabatic effective liquid water 
content (continuous curve) exc-^pt at the higher levels. The sharp cutoff 
at the higher levels is due to the reduced collection efficiency of hail 
because of glaciation. No collection takes place above the level of the 
-^OC isotherm (10 km). 

The slope of the updraft corresponds roughly to what is known about 
the Fleming storm (cf. fig. 11 of Brownirg and Foote), The width of the 
updraft at cloud base is taken as 8 km (core width - 6 km). At a height 
of 9 km, the model updraft narrows to about half its width at the cloud 
base. 

3. CALCULATIONS AND RESULTS 

Ha i 1 s tones are initiated on embryos of diameter 0,1 cm. The embryos 
are placed inside the updraft with even spacing in the vertical and hori" 
zontal of 1 and 1/2 km respectively. The embryo placement covers the height 
range from 5 to 10 km; the latter is the height of the -40C isotherm while 
embryos of the assumed size are considered unlikely to occur below the 
former height. The growth and motion of the embryos is followed using 
analytical solutions of the equations given in Section 2. More embryos are 
subsequently placed in regions requiring greater resolution. During the 
computations embryos which move out of the computational box (fig. 2a) 
are excluded from further consideration and are designated ’’out of cloud 
particles 




3 . 1 Results for the model Fleming storm 

Some representative hailstone trajectories are shown in fig. 3. Of 
the 57 evenly spaced embryos, 38 fall out of the computational box; thus 
2/3 of the embryos are designated "out of cloud." Most of these particles 
grow only for a short time in their upward trajectory before entering the 
no-growth region above 10 km height and being expelled from the computational 
box (e.g. trajectory ?. in fig. 3). A similar feature can also be seen in 
Robertson's model (Robertson, 1975). The final diameters (at OC level) of 
these particles are less than 0.3 cm. The starting region of these embryos 
is indicated by the shading in Fig. ^4. 

Ten of the 19 embryos, which fall inside the computational box, dis- 
play a simple arcing up and down trajectory (e.g, trajectory I, fig. 3). 

These particles start to the left of the shaded "out of cloud" region 
(see fig. 4). None of these particles attained a diameter greater than 
0.5 cm at the OC level. These particles should melt completely before 
reaching the ground. The growth time of these particles is about 10 to 
20 minutes (shown by the full curves in fig. ^) . 

Hail of diameter about I cm and larger originates on embryos which 
start in a narrow region on the right side of the updraft (trajectories 
3 and k, fig. 3). This region is small; it contains only 9 of the 57 
embryos. Most of the embryos starting in this region execute one or more 
recirculations. The number of recirculations increases as the embryo 
starting position approaches the right edge of the updraft, and the particles 
with the largest number of recirculations start at the 9 km height, which 
is the height at which the updraft turns over. Many of the particles spend 
over half of their time in the recirculation region and grow to half their 




Figure 3 . Representative particle trajectories in the model 
Fleming atorm^ The diameter and elapsed time 
(both at the OC level) are indicated. 
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Figure 4. (Top) Contours (dashed) of hail diameter, am, and elccpsed time, 
minutes (solid). (Bottom) Fall-out location and diameter of 
hailstone. The diameter and time refer to the OC level. 



final size there. The final diameters of these particles range between 
1.7 and 5.7 cm, and the growth times range between 40 and 60 minutes (see 
fig. k) . The maximum size (5.7 cm diameter) is in good agreement with 
the baseball size hailstones reported for the Fleming storm during its 
quasi-steady phase lasting over 3 hrs. 

Histograms of final diameters at the OC level and elapsed times 
(fig. 5 ) show the number of embryos producing the hail of indicated size 
and the time taken for the hail to reach that level. The elapsed time 
histogram shows a bi modal distribution; the low mode is composed of embryos 
originating on the left side of the updraft while the high mode is due to 
the particles originating on the right side of the updraft. The final 
diameter histogram is similar to the elapsed time histogram. A large fraction 
of the embryos produce hail of diameter less than I cm. These embryos start 
on the left side of the updraft. There is a drop in particle number around 
1 cm diameter and then a lower peak occurs. This peak is due to embryos 
starting on the right side of the updraft. 

An interesting correlation was found between the fall-out location and 
the final diameter of the hail. The major portion of the hailfall begins 
about 2 km behind the updraft with exclusively large hail. The largest 
hail size decreases with increasing distance behind the updraft edge (see bottom 
portion of fig. A). Figure A (bottom) also shows a tendency for hail of diameter 
greater than 1 cm to bunch together into groups behind the updraft. The bunch- 
ing persisted even when the embryos were distributed with a finer spacing. It 
is possible that the bunching is an intrinsic property of particle sorclng 
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(a) Relative numbers of embryos (initially j uniformly distributed 
over the updraft^ see text) csid their final diameters. (b) 
relative numbers of embryos and elapsed time. The diameters and 
times refer to the OC level. 




in patterns of sloping updrafts. With a moving storm, the bunching could 
perhaps give rise to a streaky pattern of hail fall at the ground. Towery 
and Morgan (1977) and others have reported streaky patterns of hall at 
the ground. 

3. 2 Variations on the model Fleming storm 

To test the sensitivity of the results to the mode! input we 
repeated the calculations changing the values of some of the parameters. 

In the results presented in figs. 6 and 7, the model is similar to the 
model Fleming storm, except that the maximum vertical air velocity has been 
reduced from AO m s to 30 m s and 20 m s respectively. 

The case with the maximum vertical velocity of 30 m s ^ is similar 
in many respects to the model Fleming case. However, the region of the 
updraft contributing to "out-of-cloud" particles decreases, fhe maximum 
hailstone size does not change significantly. The case with the 20 m s * 
maximum vertical air velocity has a much smaller "out-of-cloud" region. 

This case has a recirculation-mode region which gives rise to large hail- 
stones as in the case of th.j model Fleming storm, But it also shows another 
region from which large hailstones originate. This region is on the right 
side of the updraft near its base and yields hailstones with maximum 
diameter of about 4 cm in about 25 minutes. The trajectories of embryos 
starting in this region is an up, across and down arc with a small loop. 

The hailstones are able to balance in the high liquid water content region 
and grow rapidly. This case also yields 1 cm diameter hailstones in 15 
minutes and 2 cm diameter stones in 20 minutes via arcing trajectories. 

The results for this case bear a great deal of resemblance to the results 
of Paluch (1978). 
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Figure 6. Simitar to Fig. 4 but with the maximum vertical air 
velocity reduced to 30 m 8 





Figure 7. Sinvitav to Fig. 4 but with the maximum vertical air velocity reduced to 





A calculation was also performed by changing the width of the updraft 
of the model Fleming storm from 8 km to 6 km (the width of the transition 
region was kept the same). In this case, the elapsed time for trajectories 
was similar to that in the case of the model Fleming storm. However, the 
maximum hailstone diameter decreased from 5.7 to 3.6 cm. This decrease 
was due to a decrease in the number of recirculations of hailstones. 

The general results of this case are similar to those of English (1973). 

3 . 3 Particle trajectories 

Three types of particle trajectories have been encountered In the 
model Fleming storm depending upon the starting position of the embryo. 

These trajectories are summarized In fjg. 8e and compared with the results 
of Paluch ( 1978 ), fig. 0(a) through 8(d). Paluch used three-dimensional 
motion fields obtained from Doppler radar observations and calculated particle 
growth and trajectories assuming the motion field to be steady; the micro- 
physical parameters were also prescribed and assumed to be steady. The 
trajectory types obtained by Paluch (fig. 8a, b, and c) are basically 
similar to the trajectory types 1, 3 and k (fig. 8e) obtained by us. Hence, 
it appears that our conclusions, for the model Fleming storm, might be 
applicable more generally to storms with sloping updrafts, 

k. SUMMARY AND CONCLUSIONS 

Analytical solutions to the equations of growth and motion of hail- 
stones have been obtained for steady-state linear variations of the vertical 
air velocity and effective cloud water content (product of collection 
efficiency and cloud water content) with height. Fairly general models 
of quasi-stcady storms with sloping updraft can be constructed by piecing 
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• Figure 8. A swmsry of particle trajectories obtained by 

Falucht 1978 (a, bj c and d) , and those obtained 
here (e), llote that d refers to a plane perpen- 
dicular to that of a, b and c. 



together sections with piece-wise linear variations of updraft and effective 
cloud water content. In this way a model of the Fleming storm has been set 
up and the growth and motion of hail embryos in the storm computed. Three 
distinct situations arise depending upon the starting position of the 
hail embryo. Embryos starting in the middle portion of the updraft 
travel far from the main updraft and attain only small sizes; most of these 
hailstones are not expected to survive the fall from the OC isotherm to 
the ground. Embryos starting on the upshear edge of the updraft undergo 
an arcing up, ovur and down trajectory and give rise to hailstones of 
diameter less than 0.5 cm; these particles also are not e.xpected to survive 
the fall from the OC level. Embryos starting on the downshear edge of the 
updraft give rise to hall of diameter greater than 1 cm, the maximum 
diameter being 5.7 cm. These embryos execute recirculations, with the 
number of recirculations and the final hailstone size increasing with the 
proximity of the embryo starting position to the turn-over location of 
the updraft. The existence and location of this region (temperature, verti- 
cal air velocity) are believed to be important for the growth of large hail. 
It is believed that the three types of particle trajectories identified here, 
and the associated growth patterns would obtain in most cases of quasi- 
steady sloping updrafts. 
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ABSTRACT; 


A simple model of the evolution of particle size distributions by 
coalescence and spontaneous and binary disintegrations is formulated. 
Spontaneous disintegration involves single particles, while coalescence 
and binary disintegrations involve pairs of particles. Analytical 
solutions for the mean mass of the distribution and the equilibrium 
size distribution are obtained for the case of constant collection 
kernel and disintegration parameters. At equilibrium, the forms of 
the size distributions are identical under the action of coalescence 
and either or both disintegration processes; the particle concentration 
is proportional to the total mass concentration (M) and the mean mass 
of the distribution is independent of M when only coalescence and binary 
disintegrations are operative. At small values of M, the effects of 
spontaneous disintegrations dominate over those of binary disintegration'' 
while the reverse is the case at large values of M. Some of the findings 
of the present simple model are in qualitative agreement with the results 
of numerical calculations of the evolution of raindrop size spectra with 
realistic formulations of drop coalescence and breakup. 



I . INTRODUCTION 



Coalescence and breakup are among the important processes gr -erning 
the evolution of precipitation particle size spectra. Coalescence alone 
tends to increase the concentration of the larger particles^ and decrease 
that of smaller particles. Particle breakup opposes this tendency of 
coalescence and it is conceivable that a combination of coalescence and 
breakup would produce stationary particle size distributions. Indeed 
this was found to be the case in a number of numerical studies of the 
evolution of particle size spectra. Srivastava ( 1971 ) found equilibrium 
size spectra resulting from the action of coalescence and spontaneous 
breakup of raindrops. Studies of coalescence and spontaneous and 
collisional breakup (e.g. Gillespie and List 1978 , Srivastava 1978 and 
Young 1975 ) also showed the development of an equilibrium size distribu- 
tion for raindrops. Furthermore, it was found that: (i) the effect of 

collisional breakup dominates over that of spontaneous breakup at large 
precipitation mass contents (M) , while the reverse is the case at small 
M, and (ii) the equilibrium size distributions resulting from coalescence 
and collisional breakup are parallel to each other, or more explicitly, 
at equilibrium, the particles concentrations are proportional to M. 

In this paper, we set up a simple formulation of particle coalescence 
and breakup which admits of analytical solutions. These solutions shov*^ 
explicitly the approach to equilibrium and the relative importance of 
spontaneous and collisional breakup in determining the particle size 
distributions. The solutions also reproduce qualitatively some of the 
results found in the numerical calculations cited above. 



2. MODEL 6 EQUATIONS 


We consider a discrete size distribution, the particle masses beinci 
taken as the integers. The process of coalescence cannot produce particles 
of non-integral mass. The process of particle breakup will be formulated 
so that It too cannot produce particles of non- integral mass. 

The usual formulation of particle coalescence will be adopted, that is, 
the rate at which particles of masses j and k combine to increase the concen- 
tration of particles of mass (j + k) will be taken as K(j,k)p.p. , where p. is 

J J 

the concentration of particles of mass j and K(j ,k) is the collection kernel. 

Two modes of particle breakup will be considered. First, spontaneous 
breakup of particles of mass k will be assumed to occur at a rate where 

a, the spontaneous disintegration parameter, can in general be a function of 
k. The second mode of particle breakup, involving a pair of particles, is 
similar to collisional breakup and will be referred to as binary breakup. 
Binary breakup involving particles of masses j and k will be assumed to 
proceed at a rate BPjPj^ where the binary disintegration parameTer 3 can in 
general be a (symmetrical) function of j and k. All fragments resulting from 
disintegrations will be assumed to be of unit size. 

The time rate of change of particle concentration is given by 


d-r 'I jSi Pj 


(I 


- P|^ jl, K(j,k)pj - a(k)p|^-P|^ 3(j,k)pj, k > 2 

dT “ ■ jl| ^ J2 

I jll k^i (j+k)PjP|^ - Pj jlj 6(l,J)Pj 0 

For k > 2, the first two terms on the right hand side of Equation (la) repre- 
send respectively the rates of gain and loss of particles of size k by 
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coalescence, the third term represents the loss of particles by spontaneous 
disintegration and the last term represents the rate of loss of particles 
by binary disintegrations. Equation (lb) gives the rate of change of concen- 
tration of particles of unit size. The first term on the right hand side of 
the equation represents the loss of particles by coalescence. The second 
term represents the gain in particle concentration by spontaneous breakup, 
while the third and fourth terms represent the gain by binary disintegrations. 

It can be shown that Equations (l) ensure the conservation of total 
mass concentration 

N = J| l<Pk (2) 

in the system provided K and 3 are symmetric functions of their arguments. 

A general solution of equations (l) is not possible. For the coalescence 
only (a=>3*o) problem, analytical solutions were reported in the literature 
(e.g. Scott, 1968 ) for special forms of the collection kernel, namely, a 
constant kernel , a kernel proportional to the sum of particle masses, a kernel 
proportional to the product of particle masses, and a kernel which is a linear 
combination of the three foregoing kernels. In the following, we shall give 
an analytical solution of Equations (1) for a constant collection kernel K 
and constant breakup parameters a and 3* Other forms of K, and a and 3 will 
be briefly considered in Section 4. 

3. SOLUTION OF EQUATIONS AND RESULTS 
3 . 1 . Constant Ke rne I and Breakup Parameters 
We introduce the generating function 

G(z,t) = P|^ (t)z*^ (3) 

It can be shown that for constant K(=c) ,a and 3. the generating function satisfies: 
3G/at + (cN+3N+a) G = (c/2)G^+ M (a+3N)z (4) 

where 

N = N(t) P(_. (t) = G(l ,t) 


(5) 



is the total particle concentration. The following scaled variables will 


be used. 


9(z,t) 


cMt 

a/cM 

3/c 


P^(t)/M 

G(z,t)/M P,^ (t)z'' 


( 6 ) 


The mean mass of the distribution m is given by: 
m (t) = M/N(t) = l/g (1 ,t) . 


(7) 


If an equilibrium solution results, the equilibrium quantities will be dis- 
tingushed by a subscript E. In terms of the scaled quantities, Equation (k) 
may be rewritten as: 

3g/3x + (a.^ + “ 9^/2 -*■ + 3x/m)z 

We shall now consider the behavior of the mean mass and the equilibrium 
particle size distribution. 

3.1.1 Mean Hass 

Putting z=l in Eq. (8), we have for the mean mass of the distribution: 
(1/2 + 3, J + (a,,. - 3,^)m - a.,m^ 


( 8 ) 


(9) 


First, we consider certain special cases. For the case of coalescence only, 
we have; 

m (t) = m (o) + t/2 (lO) 

Thus the mean mass of the system increases linearly and the concentration tends 
to vanish with increasing time. 

If coalescence and spontaneous breakup are considered, the solution of 
Equation (9) is 

y =(yo + tanhx)/(l + y^tanhx) (ll) 
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where y and x are related to and t by Equations (15) and (l6) repsectfvely 
but with Sjf * o. The initial condition is y * y^ at x * o* It is seen from 
Equat ion (11) that as x y I ; thus as x an equilibrium mean mass is 

approached which is given by the following equation 

mg - 1/2 + (l/i» |/2a^) (12) 

As might have been anticipated^ the equilibrium mean mass decreases with 
increasing value of the scaled breakup parameter As m^ I , that 

is all particles tend to be of unit size. Equation (12) also implies that for 
fixed a/c, m^ increases with M. For H large^ or more precisely, for 

cN(o) /2a » I, nig 2 (c/2a)M 

It can be shown that the approach to equilibrium takes place monotoi/ ica 1 1 y . 

An example is shown in Fig. I. In this example mg = 215* For an initial 
m(o) » kSO , mg is approached to within 10^ at T =» AOO. Coalescence alone 
would have doubled the mean mass at T = 200. |t can also be shown that the 
equilibrium is stable, that is, perturbations of m about its equilibrium 
value tend to return m to the equilibrium value. 

If only coalescence and binary disintegration are considered, the soIu~ 
tion of (9) is: 

m(x) = m(o)l + (I + 1/23.P (1 - i (13) 

As T the following equilibrium mean mass is approached: 

mg » I + 1/23^ = 1 + c/28 (lA) 

The equilibrium mean mass tends to i as 3* “• In contrast to the previous 

case, mg is independent of the mass concentration. This suggests that at 
equilibrium, plots of the size distribution for different mass concentrations 
will be parallel to each other. This is verified later. It can also be shown 
that, as in the previous case, the approach to equilibrium is monotonic and the 
equilibrium is stable. 




F-igure 1. Tuo eacanpte-e of the variation of mean mass with scaled 
time for aoalesaenae and spontaneous and binary breakup 






In the general case of coalescence and spontaneous and binary disinte- 
grations the solution of (9) is given by (II) where y and x are related to 
m and t by 


y - (m - - i~) /L0/2a^ + - 3*) / ^ 


X. j^(l/2a,^ + 3,v/a,v) " 3^) J 


2 - 1/2 


(15) 

( 16 ) 


and y^ I s the initial value of y . As t an equi librium mean mass i s 

approached: , 

2 '/Z 

m^ » (a,^ - 3,y) /2a^ + [_(l/2a^ + 3^/a,y) + {% - 3*) 3 ^'7) 

It may be noted that (12) is a special case of (17)» obtained by putting 
3,1. = o in the latter. Also the limit of (17) as o gives equation (|i<). 

Again, it can be shown that the approach to equilibrium is monotonic and the 
equ i I i bri urn is stable . 

The behavior of is summarized in Fig. 2. The thin curves show the 
variation of m^ with a.,, for selected 3,y The thick curve shows m^ as a 
function of 3,y for a.^ = o. It is seen that m^ decreases with increases in 
both a,,, and 3^. 

It may be recalled that the scaled spontaneous disintegration parameter 

a.,. (= a^/cM) involves the mass concentration M, It is perhaps more instruc- 
tive to consider the behavior of m^ with M for fixed a/c and 6/c. This is 
done in Fig. 3. With coalescence and binary disintegration (a,^ = o) , m^ is 
independent of M (dotted curve). With coalescence and spontaneous disintegra- 
tion (3jj - o) , m^ varies approximately as M (dashed curve). The full curve 
is for coalescence and both types of disintegrations. At samll M, this curve 
approaches the spontaneous breakup curve, and at large M, it approaches the 
binary disintegration curve. Thus, at small M, the effect of spontaneous 
disintegration predominates over the effect of binary disintegrations; as M 






ilibrium meon moss 




. Ft-gitre 3. Variation of mean maas with total mass aonaentration 
for aeleoted valuea of the diaintegration parameteva. 



increases, however, the reverse is the case. A similar effect occurs 
for realistic formulations of collection and breakup of raindrops 
(Srivastava 1970). 

3.1.2 Equilibrium Particle Size Distribution 

The scaled equilibrium generating function is found by putting 
3/3t * o In Equation (8): 

g^/2 - g^ (ct^ + 3,v/m^ + 1/m^) + (a,^ + 8y</mg.)z = o 
considering the coefficients of z in (18) we find: 

’’lE “ + 3^/nig) /{a^^ + + 1/m^) 

or, using Equation (1?) 


» l/(2m^-l) 

Solving Equation (18) for g^, differentiating the result repeatedly with 
respect to z, and putting z = 1 , we find the following recursion equation 

/ P - 


( 18 ) 

(19) 

( 20 ) 


( 21 ) 


Equation (21) can also be written as: 

„ _ (2k-2)l 1 ^"’e ■ "’e^ 

^kE ■ k! (k-l) ! „2k-l 

^ (m^- 1/2) 


k-l 


2k- 1 


( 22 ) 


From Equation (22), it is seen that the equilibrium size distribution depends 
explicitly only on m^. The dependence on the parameters and B,y is implicit 
through the dependence of m^ on those parameters. In other words, the form 
of the equilibrium size distribution is the same for (!) coalescence and 
spontaneous disintegration, (ii) coalescence and binary disintegration, 

and (iii) coalescence, spontaneous and binary disintegration. Consideration 
of the particular forms of m^ shows further that in case (ii), the equilibrium 
distributions for given 8,^ and different M are parallel to each other, the 


particle concentrations being proportional to M. This is no* so in cases 

(i) and (iii). This result is understandable. |n the case of coalescence 

and binary disintegration, all terms in the equation for the equilibrium 

distribution are proportional to the products of particle concentrations. 

Therefore if one solution pj^^ for a mass concentration M has been found, 

A is a possible solution for a mass concentration AM where A is a 

constant. In the case of coalescence and spontaneous breakup, however, 

the particle concentrations enter Into the equation both linearly and 

quadrat i cal ly and hence we do not anticipate size distributions which are 

parallel to each other. Gillespie and List (1975) and Srivastava (1978) 

found similar behavior with realistic formulations of coalescence and 

spontaneous and collisional breakup. 

For large k, we can approximate in Equation (22) by 

(m, - 1/2) 

P = __L_ 

kE , 2 \ 

(m^ - m^) 

Lushnikov and Piskunov (1977) considered the problem of coalescence and 
spontaneous disintegration and obtained equations similar to our Equations 
(21), (22) and (23). Here we have shown that these results are also valid 
for coalescence and binary disintegration and a combination of coalescence, 
spontaneous and binary disintegration processes provided the appropriate 
m^ is used in Equations (20) through (23). 

Fig. i* shows a plot of the equilibrium distribution for m^ = 50, 
computed from Equation (22). Although indicated as a continuous curve, 
the distribution is discrete. On this figure, the asymptotic form (23) 
is indistinguishable from the exact solution for k ^ 8. For cloud physical 
purposes, it is perhaps more interesting to see the distributions plotted 
as log (p|^. k against k ^ (Fig. 5) which corresponds essentially 


m^ - m^ 


(m^ - 1/2) ‘ 




/i 


scaled particle concentration, 



Figure 4. Scaled equilibrium size distribution (full line) 

for an equilibrium mean mass = 30 compared with 
the asymptotic solution (dashed line). 



to a plot of log N(D) vs D where D is the particle diameter and N(D) the 
particle concentration density. In Fig. 5, the scaling has been removed 
and the equilibrium distributions for coalescence and spontaneous breakup 
(a/c =0.1 B = 0) have veen shown for the indicated mass concentrations. 
Curves identical in fo rm to those in Fig. 5 would be obtained for other 
cases provided the m^ are identical. As an example, let us consider the 
M = 2 curve in Fig. 4. This curve has m^ = 3-70 (Eq. 12). Under the 
action of coalescence and binary breakup, the same would be obtained 
for 3/c = 0.185 (Eq. 14). Hence the M = 2 curve in Fig. 5 is also the 
equilibrium distribution for coalescence and collisional breakup for 
B/c = 0.185, and M = 2. However, if M is changed with 3/c kept fixed, 
the equilibrium distribution will be displaced parallel to itself. 

2 Other forms of the co 1 1 ec t ion kernel and breakup pa rameters . 

Solutions have also been obtained for constant breakup parameters and 
(a) collection kernel proportional to the sum of particle masses and (b) 
col'fcct'on kernel proportional to the product of particle masses. A 
detailed discussion of these solutions will not be given here. Only certain 
new features which emerge in the behavior of the particle mean mass are 
summarized below. 

As the power of the collection kernel increases, an equilibrium mean 
mass is not always possible. In case (a), an equilibrium mean mass is 
always obtained with coalescence and spontaneous breakup. With coalescence 
and binary breakup, however, an equilibrium M occurs only for 3,.,>1 and the 
equilibrium value is independent of the mass content M. For 3 ,l< 1> coale- 
scence dominates, that i s m as x 

In cases (b) , m exhibits three distinct kinds of behavior. In a region 
of relatively large and 3y, an equilibrium m is always attained. In a 
region of small and 3,.., coalescence dominates in a finite time x^. 




breakup . 





i .e. , m as T -+ T . In a third region of sma| I a. but relatively large 
3,.,, the behavior of m depends upon the initial condition. For small initial 
m an equilibrium m is attai ned wh i le for large initial m, coalescence 
dominates. Moreover^ two equilibrium values of m are possible. The smaller 
m is stable with respect to small perturbations, while the equilibrium point 
with the larger m is unstable. In the latter case, a small positive pertur- 
bation causes coalescence to dominate while a small negative perturbation 
causes the m to move to the stable equilibrium point. In this case also the 
equilibrium m is independent of M for coalescence and binary breakup. 

Besides the above, it is also possible to obtain solutions for certain 
forms of variable breakup parameters. For example, we can consider 
«(j ,k) = cjk 
a(k) = a^k 
B(j,k) = e^jk 

where c, a and 3 are constants. A reference to Equations (I) will show 
o o 

that in this case, at equilibrium, the quantity jp^ obeys the same equation 
as the quantity pj in the case of a constant kernel and breakup parameters. 
Hence the solution to this problem is immediate. 
k. CONCLUSIONS 

The equations for a simple model of particle coalescence and breakup 
have been formulated. Analytical solutions of the equations have been 
obtained for a constant collection kerne) and constant spontaneous and binary 
disintegration parameters. The behavior of the solutions is qualitatively 
similar to that of numerical solutions of the equations for the evolution 
of raindrop size distributions for realistic fcrmulations of coalescence 
and breakup. An equi I i br i urn part i cle size distribution is attained. For 
given breakup parameters, the effect of spontaneous disintegration dominates 
over that of binary disintegration at small mass contents, M, while the 



reverse is the case at large M. With coalescence and only spontaneous 
disintegration, the equilibrium mean mass ra^ increases with M while it 
is independent of M when coalescence and only binary disintegration are 
considered. In the latter case, the equilibrium particle size distribu- 
tions for different M are parallel to each other. 
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DOPPLER RADAR STUDY OF A REGION OF WIDESPREAD 
PRECIPITATION TRAILING A MID-LATITUDE SQUALL LINE 


Ttnw J« MftMjkA tad 1— nih C* Srlvascava 
D«p«rOHac ot GaoplipsicaA Seimcmm 
tbm Oaivanicy of Chlcafo 
<aiicoin» IlZiaola 


Loof Xlvad incaoM oquAll lln— mwm 
wuitalnad nd pcopasacad by elm loltloeloo of 
oov eoanroetlon on choir looiilnt facoo. Am eho 
chuadorocopM oMoelot^i vich cho om eoovoecioo 
(lACttro and diooipaco, ch«y aovo co tho rur of 
eho propagotiag lino. Undor oortaln coodiCloao 
not wall undorocood mt cho prooonc c1jm» cho 
i&vlXo and dobrlo aaoociacod vlch cho dlooipocing 
chundoroconM aid aorgo cogochor co prtKlttco a 
largo ■oaoacala roglon (aroa of doop 

and rmlaci^iy boaogonaoua proclpitacion vhi.ch 
craiXa tha llao of arclvo eonoaccioa. Such 
aqnall Uaaa have boon daacrlbod by Zipaar (1969, 
1977) and Bouca (1977). 


Fig. la ia a PfX dlaplay of radar 
roflacciTlty factor at ZZI5 and ohowi tha aquall 
lina cc a ataga whan, tha vidaspraad praeipicatloo 
to cha caar of tm aetiwa lina waa juac twgijmlng 
CO baeona astanalwa^ By 2346, cha eoowMtiva lina 
had wa a iwn a ri con a i d ar ably , and cha araa of vida> 
apraad pracipitation had aniargad, axtanding aloft 
behind cha lina for a diatanca of 140 ka (Fig. lb). 
By chia ciaa alao cha lina of Incanaa activicy bad 
tweona datachad fron chrA araa of inconaa pracipi- 
cation by a gap of oaar 20 kn (Fig. lb). Slnca 
light, vidaapraad pracipitation conclnuad for two 
■ora hoora, cha aaciaconca of vaaMcaia aacanc co 
prowida addicional condanaata aaana likaly. 

3. BZfBSDED HKXBOD OF FAD ABALTSIS 


Obaanraciana and chaoratical cooai- 
darationa auggaat chat cha region of cralliog 
pracipitation nay ba dynanlcally aetiva and 
contain aaaoaeala rartical air aotiona. Thia la 
auggaacad by cha long Ufa of cJha ragion and eba 
fact chat Ic aoaaciaaa accoimca for a largo 
fraction ("O.A) of chci total praciplcaclon of 
cha aquall lina syacan (iiousa, 1977). Phyal- 
cally, ic can ba saan why vortical aotiona occur 
ia cha ragion of vidaapraad pracipitation if 
oondiciona ara faworabla for avaporation la cha 
lowar layara; cha cooling aaaociatad vlch cha 
araporation (iwrhapa anguinuHi by cha aaltlag 
of ica) can glva rlaa Co praaaura falla and 
convargaoca la cha aiddla cropoaphara, which 
should raault in aacanc la cha upper lavala and 
daacanc in cha lowar lavala. Tha aacanc could 
load co cha production of aora pracipltacloa and 
parhapa a aalf-auatalaing ayacaa. 

In chia paper, vi praaanc Doppler 
radar obaarvationa of a ragion of vidaapraad 
pracipitation trailing aa Intanaa aquall Una. 

The obaarvationa w a re obtained during Project 
NIMROD, conduecad la lUinola during Hay and 
Jiaia 1978. Throe Doppler radara (tha NC4£ CP3 and 
CPA and tha Uolvaralty of Chicago * lUlnola 
Staca Wator Survey CHILL) woro uaad baaldaa 
surfaca laatrunantatlon (tha NCAR PAM ayacan and 
cha ISIfS rainganga natwork) and fraqwanc rawln^ 
aonda aacanca. Bara va ahall praaanc obaanra* 
clcna from only one radar (CF3) , which wart^ uaad 
CO calculate hlgb-ranoludon vertical profii«n 
of partlcla fall a^ad, borlaontal dlvorganca, 
and vortical velocity by a andaat ontanalon of 
Cha FAD nacbod (BrowiOns and VIoalar, 1966). 

2. BEIEF DESCKIPnOH OF THE SQUALL LIME 

Tha band of Intanaa praciplcaclon 
froa cha aquall lina, followed by cho rogioo 
of light, vldoaprood proeiplcatlon, paaaod through 
northom Ulinola boewoan 2200 CDT (all cinoa ara 
Cancral Daylight Tlaa) , 17 Juna and 0200, 18 June 
1978. 


In a ragion of horlaoncally atratlfiad 
pracipitation wlch linear varlatlona of cha hori- 
cental wind, the Dopplar radial voloclcy at a 
eonacanc alavt^clon angle a, and eonacanc slant 
range R (horlaoncal range r) ia ralacad co tha 
aaiouch angle ^ (nsaaurad clockwlaa fron north) 
by cha aquation (Browninf and Wanlar, 1968) 

F • a^ -f a^coay ^ a^ainA b^coalA b^ainl^ (I) 


where 


vcoo a 


ucoe a 


( 2 ) 

(3) 


Bore Div la cha average diverganca of cha horlaon- 
cal wind over che area of cha FAD clrcla, 7 la cha 
average particle fall af»ad (poaiclva downwarda) 
over Che clrcusafarence of cha FAD clrcla, and u 
and V ara raapactlvaly cha aaaewerd and northward 
conponanca of cha wind velocity over cha radar. 

Tha coafflcianca b. and b2 in (1) are relaead co 
cha dafomatlon of^cha horiaontal wind and will 
HOC ba uaad hare. Tha uaual pracclca in FAD 
analyaia haa bean co caka obaervetiona of F, aa a 
function of cha asiaucb at one alavacion angle 
and CO calculaca cha coefflelanca a^,a. ,a.,b. 
either by hamonlc analyaia of cha F(<^J or by a^ 
laaac*aquarea fitting procedure. Obaervetiona in 
different range bine, correaponding co dlffnrenc 
R (and chua having different r and height h) , chan 
give cha height profllea of cha eoafflclanta. Tha 
coeffldanc a cna ba uaed co dacamlna che verti- 
cal profile of horiaontal ilvergenca, and chance 
by Incagrttlon cha profile of cha vertical air 
velocity w, provldad 7 ia known or cha alavacion 
angle a la snail enough co naglact cha laat cars 
in Equation (2) ; in none prt&vloua work, a waa 
kept sn ail so chac che latter waa cha case. 


Hare va have extended che FAD nachod, 
so chac obaarvationa of F(4)> , at a nunbar of ela- 
vatlon anglaa, can ba used co deduce boch cha 
dlvargaaca and Cha fall spaad, and ic la qoc 
oacaaaary (la fact, it la not dealrabla) co 
raacrlcc cha FAD scans co snail alavacion anglaa. 
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(2) M foiiova; 


(4) 


rb« VAD acona at diffaront oiavatloa oDtlaa cod 
b« uMd to datormliM cbo a^a ood tho quontlty 
oa cho loft hood tldo of (4) for tho rongo 
thot foil Into a aorrov tmight Intorrol. Tbo 



(cl) 



Fig, ?PT displayM of CP- 3 radar r^fisotixnty 
factor . Slack >40, hatched 30-40, dsneely a tippled 
20-30, aparaaiy stippled <20 JBZ. Cro!%e indicates 
position of radar. lAaxiimwt range is 108 km, aJ 
2225 CV7; alaration 4. 5^i squall line is to norths 
jest, (h) 2346 7DT; elevation 3,5**; squall line 
is far to southeast. In ibJ , naxima 50 km to the 
•wrth md aoutheest are associated jith the 
netting band. 


loft hood sldo of (4) LhoB tloo a dlotlnctlvo 
dopondOBco on i ohlch con bo uood to ooporeto 
tho Olv and f. booldoo roqmirlog tho onouaptlono 
■odo in tho uouoi VAD onoiTOlo (vU. , horlaontol 
unliormity of proclpltotlon and llnoor ▼oriationo 
of tho trlndo In horlaontol loyoni), tho prooont 
nothod In noro otrlngont in tomo of roquirlng 
tomporoi valfoimiry of tho phonoaonon bocouoo 
nony VAD ocono auot bo utlilood in tho onoiyolo 
proenduro. i i ow^r , tho norhod dooo contoln on 
Intomoi toot of conolotoncy through tho roquiro- 

nont of rm oxuiicit dnpondonco of (o /rcooa) on a, 

o 

In tho prooont onoiyolo, doto fron 
21 conploto oolnuthol tcono, ranging In olorotion 
ongio fron 0.5** to 84.3*, horo boon uood. Tho 
time porlod for thooo rodor acono woo fron 2343 
to 2333. Soch oalauthoi aeon hoc" opproxinotoiy 
470 oboonrotiona of Doppior miocity. Only tho 
doto for tho rongo blno within c<40 kn of tho 
rodor woro utlliiod. 

Tho doto woro first odltod to roaovo 
"ciuttor" points (obooluto roluo of rodiol voio- 
city loos thonino*^) and vniocitioo which woro 
oboorvod In tho ontonno's sido lobo duo to 
"shndoiring" of foots. Tho rodiol ▼olocitloo woro 
vmfoldod by conporison with tho bodogroph fron 
s rnwlnsondo oocont at 2336 CDT fron tho CP3 
rodor sito. A first fit of typo (1) woo nodo 
by a looot-squoroo nothod for soch oalauthoi scon 
and rongo bin subjoct to tho rongo constrolnt 
nootioQod obowo. (Approxlaotoly 2300 fits woro 
nodo.) Tho Doppior volocitlos woro thon "cloonod" 
by ronovlng volocltlss which dovlotod fron tho 
fit by noro than tvlcs tho standard error of tho 
first fit. Tho "ciaoning" hod tho offset of 
rono^rlng spurious obsonrotions rsthor chon 
nstoorologl colly significant data. A socond 
looot-squaros fit woo thou nodo on tho "cloonod" 
points to dotomino tho coofficionts in ( 1 ) for 
soch sloTOtion sngls and rongo bin. 

Tho aoat analysis stop conslstsd of 
stratifying tho s^s In narrow height intorvolo 
(300n). A loaot squoroo fit, according to Eq. 

(4) , woo porfomsd to dotomino Div and F for 
ooch of the loyoro. An Integration was thon 
errriod out to dotomino tho vortical air volocicy. 
Tho dotoils of this aro doscribod lotor. 

For tho soko of brovity, wo hovo 
onittod, iai tho above, certain details of tho 
analysis procoduroo. Thooo are: (i) height and 

olovotioQ ongls corroction for ray-curvoturo (tho 
latter corroction woo a snail fraction of a dogroo 
at tho forthoot ronoo considorod) and (li) on ox- 
tonsion of tho fit of Eq. (4) to include quadratic 
height variations oi Div and F within a layer; this 
woo dons to account bettor for sharp variations, 
such os through tho oolting layer. 

4. DISCUSSION OP RESULTS 

4.1 Winds, tonporsturo, and humidity 

Vortical profiles of wind speed and 
direction obtained from tho VAD analysis In tho 
region of wldosprood precipitation behind tho squall 
lino and from a rowinsondo launched fron tho CP 3 
radar sito at 2336 aro shown In Fig. 2. Vortical 
profiles of tonporaturo and dov-polnt tooporaturo 
aro also shown. 
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(•) 


Fi^,2. (a) V4Ttioal profilsM of ts^^ratur^ (§olid 
tins) 2 ,id dmh^point ton^ratur^ (dashsd tins) from 
rauinMo^.ds §omding» (b) vortioal profits^ of 
'jind dirtation (§olid lino) and lAnd epo^d (doMh^d 
lin§) from ths VAD (malyHs; also from rajinsonds 
soundutg ( ^ md z). Ths vsrtioal axis is alti^ 
tJ4ds abovs ssa Isvsi, Tromd Isvsl uas 226 m. 


Fig, 3, VsTxioal profilss of hydromstsor fall^ 
spssd from ths VAD analysis (solid lifts) ca%d from 
radar data oollsotsd in osrtioally •pointing mods 
(VPMJ (^) and of radar rsflsotivity faator from 
VPH (dashsd I ins). 


Tha ctMrsAl ttructura of ths stmos^ 
ptMro rsMabiad ch« soimdliifi b«hlnd troplcAi 
liiMt da«crlb«d by Zipssx (1977) and bnhind an 
Oklahona aquail Una daacrlbad by 0|ura and 
Liou (1980). Balov 4 ka, cha aounding vaa 
cooalatanc vlch cha occurranca of avaporatlon 
of praclplcacioo and onaacuracad daacanc (Laary, 
1980). Abova 4.0 kn, cha aCaoaphara vaa aatura- 
cad vlth raapacc co lea. Baevaan 4.0 and 8.7 ^n 
and abova II. 0 kn, cha acanaphara vaa acabla. 
Baevaan 8.7 and II. 0 ka« It vaa coodlciooally 
laacabla. Tha layar froa cha aurfaca co about 
I ka vaa vary acabla. Tha vac bulb potanclal 
caaa>aracura In chla layar vaa 292 E. In cha la- 
yvT froa I co 4 kn Ic vaa noacly baevaan 290.3 
and 291.5. Thla concraac aupporta Zipaar'a (1977) 
conclualoo chac cha air In chaaa evo layara haa 
diffaranc orlglna. 

4.2 Hydroaacaor fail apaad and raflac- 

civlcy factor 

Tha vartical prof 11a of hydroaacaor 
fall apaad obcalnad froa our axtanaioo of cha VAO 
aachod la shovn in Fig. 3. Tha profllaa of avaraga 
fall apaad and raflacclTley factor aaaaurad vhlla 
cha radar vaa in cha varclcally pointing aoda 
during cha 5.5 aln iaaadlacaly folloving cha VAD 
acana ara alao ahovn. Tha evo tact of fall apaada 
ara in good agraaaanc; aapaclally notavorthy la 
cha agraaaanc through cha aharp cranalcloo In cha 
laalclng zona baevaan cha alclcudaa of about 3 and 
4 Ka. Thla agraanaoc aupporta our aachod of aapa- 


racing cha fall apaad and dlvarganca froa cha 
"avaraga" of cha VAD valoclclaa (a^) and alao 
glvaa cradanca co cha aaauapciona undar lying ch«» 
uaa of cha aachod. 

Froa II km dovn Co 6.5 ka, cha fall 
■paada ganarally Incraaaad froa I co 2 aa* . 

Thaaa fall spaeda ara typical of anov. Tha In- 
craaaa of cha fall apaad and of cha raflacclvlcy 
factor laply cha occurranca of partlcla grovch 
by aoaa or all of cha folloving procaaaaa: 
dapoaltloo, riming and aggragaclon. Baevaan 6.5 
and 4.2 in, tha Incraaaa of fall apaad vaa 
arraacad, and cha rata of Incraaaa of raflacclvlcy 
factor alao dlmlnlahad. Tha formar could hava 
baan dua co cha format ion of larga aggragacaa, 
vhoaa fall apaada ara a vary vaak function of 
chair mass, Tha latcar could parhapa hava baan 
dua co cha caanaclon of aoma of cha grovch 
procaaaaa oparaclva abova 6.5 km. 

Tha fall^spaada incraaaad rapidly 
from about 2 to 9 aa~^ Ir cha malting ragloo and 
chan dacraaaad dovnvarda co approxlmataly 
7.5-8 aa*^. Tha raflacclvlcy factor Incraaaad 
sharply from about 28 to 41 dBZ and than fall 
again to 29 dBZ. Thla Impllaa a prapondaranca of 
conaldarabla partlcla aggragaclon abova, and braak- 
up balov, tha paak of cha raflacclvlcy curva. Tha 
paak fall spaad of about 9 aa*^ vaa unusually high 
for a iialclng-band situation. It Impllaa cha 
aalacanca of rathar naaalva parclr.Iaa In cha forma 
of aggragacaa or rlmad parclclaa; chaaa parciclaa 
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OlVOaiCS 


Fig. : ^^rtiaaL profits of horisontal divsr^ 

gsnos from VAD (malysis. dashsd tins 

rsprsssnts jn e x trapolation of ths profits to 
oloud top. 


$ 


could «vm hsy bMo d«rlv«d fro* ch« d«brl« of 
th« lataaM convict loo. Tho docroaoo of cho foil 
voioclty CO « 7. 3-d in lowor Uvoio U coo- 

•loconc vlth coociDuod partlclo brookup «nd 


raindrop slzo disc rlbuc loon having ch« Monurod 


rtflacclvltr faccort. 


4.3 Horizontal dlvarganca and vortical 

air valoclty 


Tho vortical profllo of tho horizon- 
tal dlvorgonco obtalnod froa tho VAD analyola la 
thoim In Fig. 4. A diocroco fora of tho 
approalaoto continuity aquation 

Dlv 3 (j3w)/3h • 0 

vaa Intagratad to obtain tho vortical air volo- 
clty w. Tho air donalty o wao takon fro* tho 
rovlnaondo aacont. Throo altomatlvo boundary 
coodltlooa gava tho throo vortical air voloclty 
proflloa ahoon In Fig. 3. For curvo (A), w wao 
takon to bo 0 at ground lovol. Thla wao juotl- 
flod by tho oboorvod loir-lovol wlndo and tho 
known torraln of tho aroa. Curvo (B) woo ob- 
calnod by aaaualng that v • 0 at h • Ll.l ka, 
tho top of tho radar data uood to obtain tho 
dlvorgonco profllo. Abovo thla holght, tho 
radar data voro apotty, and It wao difficult to 
calculaca a rallablo va3ua of tho horizontal 
dlvorgonca. Howovor, tho cloud did aatand ^oova 
thla lovol. Thla la aupportod by tho proaonco 
of upward notion at thla holght in curvo (A) and 
by satolllta data. Tbo cloud top wao aatlnatod 



fig. S. Vertioal profitss of 'Jsrtioal zir 
osloaitiss, oorputed from intsgrations of ths 
divsrgsnos profits 'Fig. 4). 


to bo at h • 12.7 k* fro* Infrarod aotolllto 
Inagory. Tho third Intagratlon thoroforo 
aaauaod w ■ 0 at h ■ 12.7 k* and producod curvo 
(C) of Fig. 3. For thla intagratlon, tho dlvor- 
gonco profllo of Fig. 4 waa oztrapolatod ilnaorly 
to hlghorlovola, giving Dlv - 3.1al0*^a*^ at 
h • 12.7 k*. (Thla aztrapolatlon la shown by tho 
daahod lino in Fig. 4) At thla tl*o, wo bollovo 
profllo (C) to bo noaror tho actual vortical 
notion than althor of tho othor proflloa. 

It nay bo ro*arkod boro that snail changoa in the 
aaouBod boundary condition at a low lovol glvo 
largo changoa in tho intagratad w at high lovola, 
whllo saall changoa in tho aaauBod boundary con* 
dltlon at a high lovol glvo rlao to only taall 
changoa in tho co*pucod w at low lovola (Bohno 
and Srlvaatava, 1976). 

-4 -1 

Tho nazliM convorgonco , -1.3x10 s 
occurrod at h 3.3 ka in tho nolting zona. 

Abovo 3.3 ka, tho wind flold waa convorgont up to 
h a 9.9 ka, abovo which It bocaao rapidly dlvor- 
gont. Bolow 3.3 ka, tho convorgonco changod vary 
rapidly to a strong dlvorgont flold through tho 
nolting layor with a nax1a>ia dlvorgonco of about 
4xl0*^a*^ at 1.0 ka. In tho 270 a lavor noaroat 
tho ground, tho wind flold again bocaao ^.onvor- 
gont. Thla roault woa alao indlcatod by 
objoctlvo analyola of PAM notwork obaorvot Vona, 
which yloldod a convotgont wind fiald although 
of nuch saallor nagnltudo,-3xl0"^s**^ . 





Th« vvrtlcAl velocity profiXt* shov 
that cha cropoaphara waa characcariiad by a daap 
layar of daacaoc baXov aa4 a daap layar of aacaoc 
abova^ A >na«iu« dovnvard yaiocity of «»25 to *35 
cm occuirrad at h » 3 km. from oaar tha 0€ 
lavai (h « 3 km) to about 7 km bmighc^ chmcompu- 
cad varticaX nociocia ara amaXl (• 10 cm a*^ ) 4 uid 
parhmpa aaar tha Ximica of accuracy of tba analy- 
ala. Abova 7 km, aacant occurrad vich a 
of bacvaan XO and 33 cm a*^* Tha profllaa ara 
ccmalacaat vl>th Zipaar <X977), vha coucXudad that 
tha maaoacaXa dovadraft air pa&acrataa to vithin 
X km of, but doaa not raach, tha aurfaca. 

T'vm compuCad proflXaa of vartlcaX air 
vaXocity ara conalacamc fHrith tha profllaa of 
tamparatura and dau-point camparatura, l.a. , tha 
ragion of daacant vaa unaaturatad, vhiXa tha 
ragioo of aacant vaa aaar aaturation. 

tha profiXaa of divarganca and var* 
clcaX air valocicy art ramarkabXy conalatanc vich 
tha proflXaa of thaaa quant itiaa poatuXatad, or 
daducad from Indiracc ^.Xanca, by Zipaar (X969, 
1977) and houaa (1977) for tha ragion of vida* 
apraad pracipitation behind tropical aquaXX Xinaa. 
Ogura «id Liou (1960) anaXyt^ radioaonda data 
for an Oklahoma aquail Xina and obtainad pro- 
fiXaa of divarganca and vartlcaX air vaXocity 
in tha raar portlona of tha aquaXX Xin?i vary 
sifltLXar to thoaa obtainad hara. A aumarlcaX 
nodal (Brovr>, 1979) of tha ragira of widaapread 
pracipitation, with paramatariaatloa of tha in- 
canaa convection, aXao yiaXdad araaa of upvard 
and downward air isotiona of aiailar intanaicy 
and diatribution in haight. 

5. SUMMARY AND CONCLUDING REMARKS 

An axtauaiou of tha VAD nachod hAa 
baan davaXopad vhich anabXaa tha calculation of 
high-raaoXutioQ vartlcaX profiXaa of particLa 
fall apaad and horizontal divarganca (and, chara- 
fora, tha varticaJ air vaXocity). Tha nathod 
uciXizaa azimuthal scana of DoppXar vaXocity at a 
number of alavation angXaa. Tha taachod haa baan 
appliad to cha study of a ragioit of vldaapraad 
precipitation trailing an Intanaa squall line 
during Project NIMROD. Tha axtandad VAD annXyaia 
gave vertical profllaa of particle fall speed in 
excalXant agraaaent with cha profllaa obtainad by 
varcicalXy pointing obaarvaciona , avan through tha 
ragion of sharp changaa in tha malting layer « Tha 
compucad vertical air valocitiaa shoved general 
daacant balov tha OC laotharm, general aacant 
above h ^ 7 km, and small vertical air mociona 
in cha i!^carmadiaca layar. Tha profila of ver- 
tical air vaiucicy i.a conalatanc vich vh«c haa 
baan daducad to occur behind tropical squall linos 
and in an Oklahoma squall Ilna from indiracc evi- 
daaca and from radioaonda obaarvaciona, raapac- 
tivaly. Thla study givaa tha first direct evi- 
danca of tha vartlcaX air velocicioa in tha 
ragion of vidaapread pracipitation aaaociacad 
vich a squall linst 

Qufititinative knovladga of cha chorac- 
tariacica of tha ragion of vldaapraad pracipita- 
tion la important for a numbar of raocona. In 
this study va obtainad vertical air motiona of 
order 10 cm over tha analyzed area of radius 
40 km ( tha actual area could ba larger) . Tha 
trannportfl effaccad by such aacant ara equivalent 


to the craneporta vhich vould ba aff acted by lb 
alamanta each of radius 1 km having a varcical 
air velocity of 10 aa*^. Haaca cha aaaoacola 
vertical mociona in cha region behind cha active 
convection in cha squall lina .oanoc ba Ignored 
in conaidarntloiu of budgaca of tha squall lina. 
Tha varcical motiona could also have dynamical 
aignificoaca for tha life history of the squall 
lina. 
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